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Abstract 

We consider a model of p independent Ising spins on a dynamical planar <jr graph. 
Truncating the free energy to two terms yields an exactly solvable model that has a 
third order phase transition from a pure gravity region (7 str = — 1 /2) to a tree-like 
region (7 str = 1/2), with 7 str = 1/3 on the critical line. We are able to make an 
order of magnitude estimate of the value of p above which there exists a branched 
polymer (ie tree-like) phase in the full model, that is, p ~ 13-23, which corresponds 
to a central charge of c ~ 6-12. 


1 Introduction 

There has been much interest in studying models of conformal matter coupled to 2d 
quantum gravity and in particular in investigating the nature of the c = 1 barrier, beyond 
which KPZ theory [|]] breaks down. These models tend to exhibit a branched polymer 
phase for large enough values of the central charge c, and our aim is to shed some light 
on the extent of this phase for the multiple Ising model. 

This letter extends the work of a previous paper |2j] in which we studied a model of 
p independent Ising spins coupled to 2d gravity, the central charge of this model being 
c = p/2. The partition function of the model is 

a) 

c&g s ° 

where the sum is over some set Q of connected planar IV-vertex (/> 3 graphs. The symmetry 
factor Sq is equal to the order of the symmetry group of the graph G. Zq is the Ising 
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partition function for a fixed graph G with a single spin per vertex and coupling constant 
f3, that is, 

= S ‘ S i) ’ (2) 

Z °{S) V <«> ) 

with 

Z 0 = 2 N {cosh (3)^ . (3) 

Si is the spin on the z-th vertex {Si = ± 1 ) and the sum in the action is over nearest 
neighbours on the 0 3 graph. 

As in our previous paper we define three different versions of the model. For model I 
the set Q consists of all connected planar A-vertex 0 3 graphs. In model II the set Q is 
restricted to one-particle irreducible graphs and for model III only two-particle irreducible 
graphs are used. 

It has been proven || that in the limit of large p tree-like graphs dominate the partition 
function, for model I (see fig |] for an example of a tree graph). By solving a model in 
which the free energy of the Ising partition function was truncated to a single term we 
managed to show how this change to dominance occurs, for small (3. In this approximation 
the change was continuous and for finite p there was no phase transition to a tree-like 
region (which we would associate with the branched polymer region occurring in other 
similar models). 

The purpose of this letter is to extend the calculation by including two terms in the 
truncated free energy. We will show that for this two-term truncation of model I there 
exists a third order phase transition from a pure gravity region (with 7 ^ = — 1 / 2 ) to a 
tree-like region (for which 7 s t r = 1/2) and that 7 s t r = 1/3 on the critical line. The location 
of this critical line approximates to that in the full untruncated model and allows us to 
estimate the value of p for the point at which the tree-like, unmagnetized and magnetized 
phases meet in this model; the result is p* ~ 13-23, corresponding to a central charge of 
c ~ 6-12. In this letter we present the main details of the calculation, further information 
about the solution can be found in reference Q . 


Figure 1: Tree-like graph 
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2 Definition of the two-term truncated model 


For a given q i> 3 graph G , the Ising partition function can be written in terms of t = tanh /3 
as a high temperature expansion, 


Zg — 1 + n d l > 

i 


(4) 


where ni is the number of closed (but possibly disconnected) loops in the graph, which 
contain l links and use no link more than once. Defining pg by 


MG = lhn — logZ G , 

iV— KX) iV 


(5) 


we have for model I, 



nit + 



( 6 ) 


Defining n% to be the number of connected loops of length two (hereafter referred to as 
“2-loops”) then ri 2 = + \n\(n\ — 1) and we have 


MG = 


N 


nit + - ^ni^ t 2 + 0(t 3 ) 


(7) 


Suppose that we truncate mg to the first two terms of the series, denoting this new 
quantity by Mg an d consider the model with the partition function 


Zir(p) = E -e" 5 '"’. 

Gee Sg 


( 8 ) 


which we will refer to as the two-ternr truncated model. From now on Q is taken to be 
the set of model I graphs. In the limit of small (3 we might expect it to reproduce the 
behaviour of the full model. For larger (3 it will be a bad approximation to the original 
model and in particular it has no magnetized phase. However we can modify the truncated 
model slightly so that it is a reasonable approximation to the full model for large p at any 
value of f3. The key to doing this is to note that one can factor out the entire contribution 
to Zq resulting from loops of length one (“1-loops”) giving 


Z G = (l + t) ni Z' G , 


(9) 


where Z' G is the partition function on the graph G’ for which 1-loops have been cut off 
(leaving vertices of coordination number two, see fig ||). It is impossible to factorize all 
the contributions from 2-loops in the same fashion, but suppose that we consider the 
partition function 

Z" =(l + tr(l + t 2 )"S = l + ]Tn^. (10) 

i 
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Figure 2: Cutting 1-loops 







This corresponds to a model in which we count n[ the number of (possibly disconnected) 
loops of length l, but restrict each of the connected parts to be of length two or less. This 
is a good approximation to the full model if t is small or if the dominant graphs only have 
short loops in them (which is precisely what happens for large p — since we know that 
tree-like graphs dominate in this limit). Thus we might expect the modified two-term 
truncated model to give a reasonably accurate indication of the transition to a tree-like 
phase even for values of /3 which are not small. The corresponding / 1 q is 


Vg = n i lo g(! +t) + n c 2 log (l + r 


( 11 ) 


In this letter we will solve exactly this modified truncated model, defined by equations 
(|sj) and @ , and relate the results to the untruncated model. For technical reasons it 
is much easier to solve the model if we use 0 3 graphs which have a “root”; that is, in 
the terminology of A0 3 theory, we are going to use 1-legged Green functions rather than 
vacuum diagrams - this will not change the free energy and should not affect any of our 
conclusions. In this case we have 


Zn(p ) 

= £ expp m log(l + t) + 112 log (1 +1 2 ) 

Geg 

OO OO 

(12) 


= £ £ g^(N, ni ,n c 2 ) v n ^y nC 2, 

ni=0 ri 2=0 

(13) 

where we have defined 


v = (i + ty, y =(i + t 2 Y 

(14) 


and Qr^ (IV, n \, u^) is the number of IV-vertex rooted model I graphs with ri\ 1-loops 
and ri 2 2-loops; the subscript “r” indicates that we are using rooted graphs. The grand 
canonical partition function is 


OO OO OO OO 

Zr= £ e~» N Z N (p) = £ £ £ gP(N, ni ,r^) x N v n 'y n *, (15) 

N =1 N=1 m=0n 2 =0 

odd odd 

where we have put x = exp(— p). 
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3 One-term truncation of model II 


In the calculation that follows we will need the generating function for (IV, 773), the 
number of IV-vertex (unrooted) model II graphs with 2-loops. Determining this is 
equivalent to solving the one-term truncation of model II and the calculation follows that 
given in [j2j for models I and III. As in that paper one can derive a recurrence relation 
valid for N > 4, 

-271-2^ Q^ 2 \N, 77-2) + 2(?7.2 + 1) {.N, 772 + 1) = (?72 + 1) Q^ 2 \N + 2, 77§ + 1). (16) 

The generating function is defined as 

OO OO 

Z^(x,y) = E E G (2) (N,n c 2 ) x N y n + (17) 

N=4 n c =0 
even 

This is given by 

-i | OO 

Z (2 ) = -X 2 (V 3 -3y + 2) - bog/7+ E S {2) (N) x N h-‘ SN/2 , (18) 

D ^ N =4 

even 

where h = 1 — x 2 (y — 1) and Q( 2 \N) is the number of IV-vertex model II graphs, 

g( 2 \N) = E -h ~ e^ 1 °g(¥) Ar IV-i (19) 

(f)!(IV + 2)! 

4 Cayley Trees 

In the previous section we gave the generating function for n^), which is the 

number of one-particle irreducible (denoted “1PI”) graphs with a given number of vertices 
and 2-loops. If we could generalize this generating function to the case of 1PI graphs with 
777+I legs, then connecting together such graphs into Cayley trees would give model I 
graphs and allow us to keep track of both 771 and 772- 

Figure 3: (a) 1PI Green function, G'i m (b) Generating function for trees, T (c) Example 
of a rooted tree 
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The generating function for 1PI graphs with a root and m other legs will be denoted by 
G lm and drawn as in fig ^a (sometimes the root will be labelled with an “r” to distinguish 
it from the other legs). It is a sum over such graphs weighted with factors of x N y n %. Note 
that although the legs have been drawn on the exterior of the Green function some of 
them may in fact be attached to links in the interior. 

The generating function for trees will be denoted T and is drawn as in fig [jb. Fig |3|c 
gives an example of a tree contributing to T. Figure || shows that T satisfies the equation 

OO 

T= E G 'i mT m . (20) 

771=0 

It will be convenient to define a slightly different function Gi m which generates 1PI graphs 
with a root and m “exits”. An exit is defined to be a place at which a tree can be hung 
(in either of two directions) and is drawn as a link with a cross on it. For example fig |ja 
contributes to G12 and there are four corresponding contributions to G 12 . I n general 
G\ m = 2 m Gi m , but it should be noted that fig ||b can not be represented in this fashion 
and must be excluded from G12. To compensate an extra term (fig |5|c), xT 2 , must be 
added to the equation, giving 

OO 

T= Y, G lm {2T) m + xT 2 . (21) 

m=0 

This formula then correctly distinguishes between graphs which differ only in that a 
tree has been hung off a link in the opposite direction (these correspond to different 
triangulated surfaces). 

Now define a generating function for the Gi m , 

OO OO 

r (x,y,z)=J2 E E G lm (N,n c 2 ) x N y n h™, (22) 

iV>0 n,2=0 m=0 

where Gi m (N, n 2 ) is the number of 1PI graphs with one root, m exits, N vertices and n 2 
2-loops, that is, 

OO 

Gim = E E G lm (N,n c 2 ) x N y n l (23) 

N>0n,c=0 



+ 


Figure 4: Equation for T 



+ • • • 
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Figure 5: (a) Example of a graph with two exits; (b) & (c) Special case 



This is what is needed, except for the fact that the 1-loops have been ignored. The only 
contribution to ( |22| ) from 1-loops occurs in Gio and is equal to x. The 1-loops will be 
weighted with an extra factor of v and this is achieved by defining 


T(x, y , z, v ) = T(x, y, z) — x + xv. 


Hence from (21) we have that T must satisfy 

T = f (x, y, 2T, v) + xT 2 

and then the solution T(x,y,v) is equal to Z r defined in (15). 
the function T(x, y, z, v ). 


(24) 

(25) 

The next task is to find 


5 Determination of f (a;, y, 2, v) 

In this section we are going to determine T(x, y, z ) the generating function for Gi m (N, n|). 
However it is useful to find G\o(N, n^) first, as this will provide a boundary condition. 
Below we show that 

G w (N-l,n c 2 -l) = 2n c 2 G^(N,n c 2 ), (26) 

for N > 4 (N is even). 

Figure 6: (a) Rooting a model II graph (b) Special case 

(a) (b) 

€• - -C# 

Take an IV-vertex model II graph G with 2-loops. It has 2ri2 orientated 2-loops 
(ie we count each 2-loop twice, once with an arrow going in one direction and once with 
the arrow reversed). These can be grouped into equivalence classes of size s G ; each class 
consists of 2-loops which can be transformed into each other under the action of the 
symmetry group. Removing an orientated 2-loop and replacing it with a leg (as shown 
in fig |ia) produces a graph with N — 1 vertices and (n§ — 1) 2-loops. Note that we are 
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ignoring the special case fig ^Sb, for which n 2 = 3 and A n 2 = ~ 3; so the resulting formula 
will not apply for A = 2. There is a one-to-one correspondence between equivalence 
classes of orientated 2-loops and rooted graphs. The number of equivalence classes is 

£ ^ = 2»!5< 2 >(iV,„3, (27) 

G-.II S G 
(W>§) 


where the subscript on the summation means that we are summing over model II graphs 
labelled with a G , which have A vertices and n 2 2-loops. Equating this with the number 
of rooted graphs G\o(N — l,n 2 — l) gives the required result. 


Now, 

T (x,y,z = 0)=J2 
N> 0 

and substituting ( |26|) into this yields 

T (x,y,z = 0) 

where Z^ 2 \x,y) is given by ([T^). 


£ G w (N,n c 2 ) x N y n * 


n 2 = 0 


2 dzW 

x dy 


+ x, 


(28) 


(29) 


Figure 7: Adding an exit to: (a) an ordinary link (b) a 2-loop (c) a 1-loop 


(a) (b) (c) 



Having determined the number of rooted A-vertex 1PI graphs with n 2 2-loops, we 
wish to generalize this to the case in which the graphs have m “exits”. Consider adding 
an exit to a graph which has m — 1 exits and A — 1 vertices; this will produce an A-vertex 
graph with m exits. If the exit is added to an ordinary link then n 2 is unchanged (figj^a), if 
added to a 2-loop then An^ = —1 (fig [i|b) and if added to a 1-loop then An^ = +1 (fig [Tjc) . 
Ignoring the special case (fig 0c) for a moment, then we wish to make A-vertex graphs, 
with m exits and n 2 2-loops. Bearing in mind that each such graph can be produced in 
m ways (ie any of the m exits could be the one that we have just added), we have 

m Gi m (N, n 2 ) = L\G\ m —i (A—1, n 2 ) + L 2 G 1 m _i(A — 1, n 2 + 1), (30) 

where L\ is the number of ordinary links in a graph with A — 1 vertices, m — 1 exits and 
n 2 2-loops, and L 2 is the number of links lying on a 2-loop in a graph with n 2 + 1 such 
loops. One can easily show that 

3 

Li = - (A - (m + 1)) + m - 2n 2 , 

L 2 = 2(ri2 + 1). 
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(31) 

(32) 





Thus we have a recurrence relation for the coefficients of Gi m in terms of those in G 


r l m— 1 


and a formula for Gio- Note that (3C) does not apply for N = 2, m = 1 due to the special 
case (fig 0c). 

Using (|30|) one can derive the following differential equation for T(x,y,z), 


3 odT , .dF fl \ 

r si + 2i(1 - v) si - \r z+ o 

The solution is straightforward and gives T as 


^7 - ^xF + x 2 {y - 1) = 0. (33) 


f (x, y, z, v) = \/l — xz To (x(l — xz) 2 ; 1 + (y — 1)(1 — xz ) 2 ^ + x 2 z(y — 1) + x(v — 1), 

(34) 

where To(x / ,y / ) = T(x',y',z' = 0). Using (29) and (p^), this is given by 


r 0 {x',y') = x' 


x' 2 1 3 00 


+ - + - J2 0 {2) (N) x ,N Nh~ 3N/2 


h A h h 


N=4 
even 


(35) 


with h = 1 — x' 2 {y' — 1). 

6 Solution of the model 


Putting z = 2T = 2 Z r and using (^5|) gives 

xz = 2xT(x, y, z, v) + \^(xz) 2 . (36) 

As x is increased from zero there is a non-analyticity at some critical value, x c . This will 
give us the free energy of the model through y c = — logx c . In order to determine x c we 
first need to find an expression for the grand canonical partition function Z r . 

Defining e = y — 1, r/ = u — 1 and 


Y = a/1 — xz 


Q 


x 


(1 


— xz) 


3 

2 


X 

y3 > 


we can eliminate x and 2 from (^g) in favour of Y and Q giving 


(37) 

(38) 


(1 - y 2 ) = i (1 - y 2 ) 2 + 2y 3 Q [yr 0 (q, 1 + e y 4 ) + y 3 q (i - y 2 ) e + qy 3 v ] . ( 39 ) 


Now defining 


H = 




1 - eQ 2 Y A 


2 

X3 

1 — xz — ex 2 


2 

x 3 


y 2 — ex 2 ’ 


(40) 
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we have using (35) 


"o (Q, 1 + eY 4 ) = Q\ 


H 4 + H + 3H J2 <3 (2 \N)NH- 


N=4 
even 


= Q*F(H), (41) 


where we have defined a function F(H). Putting this into 

L n 4 r n 6 


gives 


YQl] F(H) + a[yQ^ (e(l -y 2 )+r/) 


o = y 4 - i + 4 


Noting that YQz = xs we can use (MG) to eliminate Y and write 


(42) 


in terms of H and 


/ = xs, 


(3e 2 ) f + |/ 4 - 4/ 3 (r? + e) - / 2 (aF(H) + + 1 = 0. (43) 

This is a sixth order polynomial in /, giving / = f(F[). 

To proceed further we require a more compact form for F(H). Brezin et al [|| give a 
formula for the 2-legged 1PI generating function T 2 (equation 64 of that paper), 

(1 - 2 t ) 2 


r 2 = 


g' 2 = r( 1 - 2 rf 


(44) 


1 — 3r 

where each vertex has a weight of g' (note that g' is equivalent to 3 g in their notation). 
Now if we take each IV-vertex model II graph, pull out an orientated link from each 
equivalence class (of size s G ) and cut it, then we get 


J2 — = 3NG^(N) 

G-.II S G 
(N) 


(45) 


2-legged Green functions. Each possible connected 2-legged Green function is generated 
exactly once, so that denoting the number of such functions by G-iiN) we have G^IV) = 
3NQ^ 2 \N). Thus defining a new symbol f, 


f = E 3NgW(N)g ,N = J2 G 2 (N)gf 


rN 


(46) 


N=4 
even 


N=4 
even 


However T 2 generates 1PI 2-legged Green functions and thus the right-hand side of the 
above equation is given by 1 — 1 — g' 2 (see fig |8|); the g' 2 subtracts off the N = 2 case. 
The critical value of g' is g' c = and the corresponding critical value of r is r c = 1/6. 
Hence, 


F(H) = H 4 + H + HT{g' = H 3/2 ) 
= HTi\g' = H 3 ' 2 ). 

We are lead to define the variable r by 

H 3 = r( 1 - 2r) 2 


(47) 

(48) 

(49) 
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(a) 


Figure 8: (a) T 2 (b) T 2 1 — 1 
(b) 



and thus 

AH 2 F(H) + 1 = (1 + 6t)( 1 — 2r). (50) 

Now defining k = f^/H/r we can rewrite ( fl3| ) as 

3e 2 k 6 r 2 + 2ek 4 r — 4 k 3 (i] + e)r(l — 2 r) — fc 2 (l + 6r)(l — 2 r) + (1 — 2r) 2 = 0. (51) 

This is a quadratic in t and can be rewritten as 

A{k)r 2 + B{k)r + C{k) = 0 (52) 

with 


m 

= 3e 2 fc 6 + 8k 3 (rj + e) + 12A: 2 + 4 

(53) 

B(k) 

= 2ek A — Ak 3 (rj + e) — 4fc 2 — 4 

(54) 

C(k ) 

= 1-fc 2 . 

(55) 


The aim of this calculation is to find the closest singularity to the origin for the grand 
canonical partition function, Z r (note that Z r = z/2) and thus we wish to express Z r in 
terms of the new variables and hence in terms of x. From (4C) we have that 

Z = — 1 — FT 2 — - 

r 2x r H\ 

and any non-analyticity in Z r is due to the singular behaviour of f /H, which is given by 

k 


(56) 


£ 

H 


(57) 


1 - 2 r 

Since the term in the square brackets is a function of x 2 it will be convenient to define a 
variable q = x 2 (bear in mind that we have already defined / = x 2//3 , so that q = / 3 = x 2 ). 
Given the definition of k, 


so that 


fc 3 = / 3 ^r = -( l-2r), 

T 2 1 


T = 


k 3 + 2 q 


and substituting this into (|^) yields 

(p )\a + 2B + AC) + (B + 4C) + C = 0. 


(58) 

(59) 

(60) 


11 









Using equations (53) to (55) this gives 


k 4 - 2 eqk 3 - k 2 


3 e 2 q 2 - 4q(r) + e) + 1 


+ 


- 4 eq 2 = 0 , 


(61) 


which is quartic in k , so that we can in principle solve for k(q). This leaves the question 
of which root should be chosen. In the limit q —» 0, we have from 
H —> 0 and r —> 0. Thus from ([57]) we require k —► 1, looking at (] 6 l]) we see that at q = 0 
it becomes k 2 (k 2 — 1 ) = 0 and the root we want is the one for which k(q = 0 ) = 1 . 

Using ( |57|) and (|59| ) we can express f /H in terms of k and q. Thus finally we have a 
formula for the grand canonical partition function, 


56) that f/H —> 1, 


Z r 


1 

2 x 


1 — ex 2 


k 


2 x r 

1 ? ’ 


(62) 


where k{x 2 ) is a calculable function given by solving (61). 


At this point the simplest course of action was to write a computer program, which 
could evaluate (^) at different values of x and, by slowly increasing x from zero, locate 
the closest singularity to the origin. The non-analyticity occurs when k(x 2 ), which is real 
for small real x, becomes complex as x is increased. 


Examination of the data for p c indicates that there is a third order phase transition 
(see fig H| which plots various third derivatives of fi c ). The location of the transition is 
most clearly seen in fig 10 , which plots the quantity 1 — 6 r evaluated at x c \ the axes used 
are t and the product pt. There are two regions one for which r c = 1/6 and the other for 
which 0 < t c < 1 / 6 . 


7 The critical line 


In this section we are going to derive an equation for the critical line in the two-term 
model and determine the exponent 7 s t r for the different regions of the phase diagram. 


It will be convenient to rewrite 


( 0 ) 


as 


k' 4 a^k 3 T 0 , 2 k 2 T &\k T a 0 — 0, 


(63) 


where we have defined the coefficients a 0 , • • •, 03 . Away from the critical line the closest 
singularity to the origin is due to a double root in this equation and thus, evaluating 
everything at x c , 

4 k 3 + 3a^k 2 + 2a,2k + ai = 0. (64) 

On the critical line we have a triple root and hence 

6 k 2 + 30 , 3 k + 02 = 0. (65) 
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Now we have three equations in four unknowns (k, g, e, rj) giving a one parameter solution. 
In general it would be difficult to solve these equations, but in this case there is a short 
cut; we know that on the critical line r c = 1/6 and so 4 q = /c 3 (from (59)). Substituting 
this into (E3) to eliminate a\ yields 


6 k 2 + lici'sk + 2 a 2 = 0 . 


( 66 ) 


Comparing this with (65) we have a 2 = 0 and 2k + 0,3 = 0. Introducing a parameter e 
gives 


e = 


d2 ’ 


k = e, 


fa + e ) “ 73 ( X + 3e2 ) ’ 


1 3 

«= 4 C 


It is easy to check that (|63|) is satisfied and to derive a formula for the critical line, 


fe e 


V = 


+ 3 - e. 


(67) 


( 68 ) 


At i = 1 , e = rj = 2 P — 1 ; the above equation gives yfe = 2(4 + \/l3) and hence the critical 
line intercepts the t = 1 axis at p ~ 7.860327. For p —> 00 , 77 ~ e 3 ^ 2 and the intercept 
with the line p = 00 occurs at the root of 


(i+t) = (i+t 2 y, 


(69) 


which is t ~ 0.6124088 (note that t = 0 is not the correct root). 

In order to calculate 7 s t r it is necessary to find out how the derivatives of Z r diverge 
as q —> q c . Putting Z = 2xZ r , 


&Z _ 2 1 [ 4 q 

dq 6 k 2 g' A: 3 . 


(70) 


where g' = 11 ; note that g 7 (g c ) = 0 and q"{q c ) / 0 in general, but that on the critical line 
q"(q c ) = 0. For the region at large p and t, r c / 1/6 and hence 


dZ 

1 

1 

(71) 

dq 

q' 

yjq - q c 

dZ 

q c )~ lstr , 

(72) 


Noting that for rooted graphs 


we have jstr = 1/2 (that is, this is a tree-like region). For the region at small p and t, 
t c = 1/6 and the first derivative of Z is finite. Differentiating (|7^) with respect to q we 

find that _ 

9^ 1 . (raj 


dq 2 yjq ~ q c 


(g - g c ) (7stl ' +1) 
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and hence 7 s t r = —1/2 (the pure gravity value). On the critical line 


dZ 1 

(q-qc)*' 


( 74 ) 


giving 7 s tr = 1/3. Thus the third order phase transition in the two-term model is a 
transition from an unmagnetized non-tree-like region to a tree-like phase. This provides 
good evidence for the existence of such a transition in the full model. Presumably as 
more terms are included in the free energy, the critical line will move somewhat and it is 
not entirely obvious whether the intercept with the line p = 00 will occur at non-zero t 
or t = 0 in the full model. 


The original aim of solving this model was to determine the extent of the tree-like 
region, in particular we are interested in finding the point at which the tree-like, mag¬ 
netized and unmagnetized phases meet (assuming that such a point exists); the value of 
p at this point is denoted p*. In previous work we made an estimate of the location of 
the tree-like region based on the region for which p c was approximately linear. It will be 
useful to compare the location of the critical line with the boundary of the approximately 
linear region (referred to as the “knee” of the solution); the linear region can be defined 
to be that for which 


dfj, c 


dipt ) 


< 5. We will take 6 = 0.0168 giving pt « 6 + 2.7 1 + 0(t 2 ) 
at the knee; this corresponds to our previous estimate that the graphs were tree-like for 
pt> 6 at small t 0 . 


In fig |lT| is plotted the line of the phase transition in the p-t plane and also the location 
of the knee. The tree-like phase is inside the linear region and the two lines are quite 
close in the area of interest (near t ~ 0.73). Using the phase line as an approximation 
to the location of the phase transition in the full model we can estimate the point at 
which it intercepts the boundary of the magnetized region. Although we do not know 
exactly where the critical line separating the magnetized and unnragnetized phases lies, 
we have that t c ~ 0.73 for the single spin Ising model ||] and we know from Monte Carlo 
simulations |ij, [7], || that t c only increases slowly with p. Taking t = 0.73 gives an estimate 
of p* ~ 23, but this is likely to be an over-estimate since the phase line changes rapidly 
with t and the value of t that we are using is too small (it corresponds to the critical 
coupling for p = 1). An estimate using our old criterion based on the knee would give a 
value of p* ~ 13. Overall then we conclude that p* ~ 13-23 and that the corresponding 
value of the central charge, above which there exists a tree-like phase in the full model, 
is c ~ 6 - 12 . 
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Figure 11: Phase diagram for the two-term model: (a) Phase line (b) Knee (5 = 0.0168) 
(c) The line t = 0.73 
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8 Conclusion 


In fig |l2] we have drawn some possible forms for the phase diagram of model I. The 
magnetized region is labelled M, the tree-like region T and the remaining unmagnetized 
non-tree-like region U. If the critical line in the two-term model were an accurate rep¬ 


resentation of that in the full model then the correct phase diagram would be fig 12 a. 
However, one expects that as more terms are included in the truncated model the location 
of the critical line will change and it is not entirely clear where point A will move to. It 
is possible that A tends towards the t = 0 axis as more terms are added in which case 
fig [j 3 b would be the correct diagram. Other possibilities, such as fig [l2| c or the absence 
of a tree-like region, can not be entirely ruled out, but seem unlikely given the results 
presented in this letter and those due to Wexler || 1C]. The two-term model has a third 


Figure 12: Possible phase diagrams for model I 





order transition with 7 s t r = | on the critical line AC and the full model may well show the 
same behaviour. It is interesting to note that the model in reference JTT| also has 7 s t r = | 
on the line separating the branched polymer and unmagnetized phases (in this model 
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7 str = \ at the point where the three phases meet). The order of magnitude estimate of 
p* ~ 13-23 (point C in fig [L2]) is similar to the largest value of p (ie p = 16) used in simu¬ 
lations so far. It is to be hoped that future simulations will manage to locate the onset of 
the branched polymer phase. One should note that our estimate of p* is not sufficiently 
accurate to exclude the possibility that p* = 2, although numerical simulations seem to 
rule this out. In conclusion then, significant progress has been made in elucidating the 
nature of the phase diagram for the multiple Ising model on dynamical planar cf> 3 graphs, 
in particular in the large p limit, but much work remains to be done before the nature of 
the c = 1 barrier is fully understood. 
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